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Abstract 

We demonstrate that both the and threshold resummations can be performed in 
the Collins-Soper resummation formalism by evaluating soft gluon emissions with infrared 
cutoffs for the longitudinal and transverse loop momenta, respectively. The reason the kx 
resummation for a parton distribution function leads to suppression in the large b region, b 
being the conjugate variable of parton transverse momentum kx, and the threshold resum- 
mation leads to enhancement in the large N limit, N being the moment of a distribution 
function, is a consequence of opposite directions of double-logarithm evolutions. The kx and 
threshold resummations for an energetic final-state jet give suppression. The switch of the 
threshold resummation from enhancement to suppression is attributed to a nonvanishing 
jet invariant mass. In the same framework we derive a unification of the kx and threshold 
resummations for a parton distribution function by requiring infrared cutoffs for both lon- 
gitudinal and transverse loop momenta. This unified resummation exhibits suppression at 
large b, similar to the kx resummation, and exhibits enhancement at small b, similar to the 
threshold resummation. 
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I. INTRODUCTION 



Recently, we have applied the Collins-Soper (CS) resummation technique [[]]] to the kx 
resummation for a parton distribution function </>(£, fey, p + ) associated with a hadron of 
momentum p^ = p + 5 fl+ [0] , which describes the probability that a parton carries longitudinal 
momentum £p + and transverse momentum kr- In this formalism the double logarithms 
\n 2 (p + b) contained in 0, b being the conjugate variable of kr, are organized. The result is a 
Sudakov suppression factor, quoted as [[J 

<f>(£,b,p + ) = exp 

where the anomalous dimension 7^ will be defined in Sec. II, and the initial condition 
</>(£) of the double- logarithm evolution coincides with the standard parton model. We have 
neglected the dependence of <ft on the renormalization (or factorization) scale /x, which is 
not relevant to our discussion here. The \x dependence denotes a single-logarithm evolution, 
which can be easily derived using renormalization-group (RG) equations. 

Equation (P has been employed to evaluate the pt spectrum of direct photon production, 
Pt being the transverse momentum of the direct photon [Q. It was observed that the double- 
logarithm exponential, after Fourier transformed to k^ space, provides the necessary k^ 
smearing effect, which resolves the discrepancy between most data sets and next-to- leading- 
order QCD {ototl) predictions [Q. This success implies the importance of the kr resummation 
for studies of hard semi-inclusive QCD processes, in which transverse degrees of freedom of 
final states are measured. 

In a series of works || the threshold resummation ||-f| for dijet, direct photon and heavy 
quark productions in kinematic end-point regions has been performed. In this formalism a 
different type of double logarithms ln 2 (l/iV), which come from the Mellin transformation 
of the logarithmic corrections ln(l — £)/(l — £) + to moment N space, is organized. They 
are important logarithms, since £ — >• 1 corresponds to the large N limit. On the contrary, 
the threshold resummation for a parton distribution function where the transverse 
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l/b p Jl/b 2 [I 



^-7Jf(««(^)) 
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degrees of freedom have been integrated out, results in enhancement: 

/"I dX 



cxp 



- I 1 dz- 

r o 1 — z J(i-z) 2 X 



7i^(a s (v / Ap + )) 



(2) 



with the initial condition <jyP> and the same anomalous dimension 7^. 

At the kinematic end points of QCD processes, energetic final-state jets may be produced, 
which also involve large double logarithms. Take deep inelastic scattering (DIS) as an 
example. The final-state quark scattered by the photon with momentum q carries the 
momentum p' = £,p + q = ((£ — x)p + ,q~,0) for the Bjorken variable x = —q 2 /(2p ■ q) = 
_q+ jp+ j n a f ram e where q possesses only the plus and minus components. As x — > 1, p' 
approaches the light cone with vanishing invariant p' 2 , and the final-state quark forms a jet. 
The threshold resummation for a jet function J(N,p'~) gives Sudakov suppression ||, 

z N ~ l -l r 1 -* dX 



J(N,p'- 



cxp 



dz- 



->y K (a s (y/Xp )) 



J(0) 



(3) 



1 — Z J(l-z) 2 X 

opposite to the enhancement in Eq. (Q). While the kx resummation for a jet function 



J{b,p' ) still gives Sudakov suppression similar to Eq. 

rv'~ dp [P 2 dfi 2 



J(b,p'- 



cxp 



7A-(a*(ju)) 



(4) 



'1/6 p Jl/b 2 fi z 

The above expression has been intensively appied to the analyses of end-point spectra of 
heavy meson decays 0. 

With the above comparision, it is worthwhile to explore the relation between the kx and 
threshold resummations, and the difference between their results for a parton distribution 
function and for a jet function. In this paper we shall show that the threshold resummation 
can also be performed in the CS framework. The double logarithms ln 2 (p + b) and ln 2 (l/A) 
are summed simply by choosing appropriate infrared cutoffs in the evaluation of soft gluon 
corrections. If transverse degrees of freedom of a parton are included, 1/6 will serve as 
a transverse infrared cutoff. Combined with the collinear logarithms lnp + , \n 2 (p + b) are 
generated. In the end-point region with x — > 1, p + /N is a longitudinal infrared cutoff, 
implying the presence of ln 2 (l/A). Therefore, in the region away from end points, we 



neglect the longitudinal cutoff, and sum \n 2 (p + b). As close to end points, we keep only the 
longitudinal cutoff, and integrate over transverse degrees of freedom of a parton in the same 
formalism. The logarithms ln 2 (l/A r ) are then summed. 

It will be shown, in the CS framework, that suppression from the kx resummation and 
enhancement from the threshold resummation for a parton distribution function is a conse- 
quence of opposite directions of double-logarithm evolutions: from the small 1/6 to the large 
p + in the former case and from p + to p + /N <C p + in the latter case. The comparision of the 
two resummations for a jet function is subtler. In the kx resummation the jet momentum 
is parametrized as p' = (0,p'~, k T ), similar to the parton momentum 0, k^), and thus 
the same suppression effect is expected. In the threshold resummation the jet momentum 
is p' = ((£ — x)p + , q~ , 0), different from the parton momentum 0, 0). The former has 

a nonvanishing invariant mass p' 2 = 2(£ — x)p • p', while the latter lies on the light cone. We 
shall demonstrate that this difference is responsible for the switch from enhancement for a 
parton distribution function to suppression for a jet fucntion. 

Once the two different resummations can be reproduced in the same formalism, it is pos- 
sible to develop a unified resummation, in which both the logarithms \n 2 (p + b) and ln 2 (l/iV) 
are organized to all orders. The unification of the k T and threshold resummations for a par- 
ton distribution function are achieved simply by retaining the longitudinal and transverse 
cutoffs simultaneously. It will be observed that the result exhibits suppression at large b, 
similar to the kx resummation, and enhancement at small b, similar to the threshold re- 
summation. The unified resummation is appropriate for studies of QCD processes, in which 
final states possess large rapidities (corresponding to large momentum fractions), and their 
transverse degrees of freedom are measured. 

In Sec. II we review the application of the CS technique to the kx resummation for 
a quark distribution function. In Sec. Ill we perform the threshold resummation for a 
quark distribution function in the same formalism, comparing each step with that in the 
kx resummation. The threshold resummation for a jet function is derived in Sec. IV. The 



4 



unification of the two resummations is proposed in Sec. V. Section VI is the conclusion. 



II. k T RESUMMATION 



Consider a quark distribution function for a hadron in the minimal subtraction scheme, 

<P(x,k T ,p + ) = J^I ^ 2 e- ixp+v ~ +ikT - yT { P \q(y-,yT)\lMm , (5) 

where 7 + is a Dirac matrix, and \p) denotes a hadron with the momentum p 11 = p + 5^ + . 
Averages over spin and color are understood. The above definition is given in the axial gauge 
n ■ A = 0, where the gauge vector n is assumed to be arbitrary with n 2 ^ 0. Though the 
definition is gauge dependent, physical observables, such as hadron structure functions and 
cross sections, are gauge invariant. It has been shown that the n dependences cancel among 
convolution factors of a factorization formula, i.e., among parton distribution functions, 
final-state jets, and nonfactorizable soft gluon exchanges, in the factorization formula for a 
DIS structure function f2|. A resummation can also be performed in the covariant gauge 
d ■ A = 0, and the result is the same as that from the axial gauge J2|. 

The essential step in the CS technique is to derive a differential equation p + d(f>/ dp + = C(p 
U, where the coefficient function C contains only single logarithms as shown below, and 
can be treated by RG methods. In the axial gauge n appears in the gluon propagator, 
(-i/l 2 )N^ v (l), with 

Nn ^ r ^±l^L + n >J±-. (6) 

Because of the scale invariance of N^ u in n, depends on p + through the ratio (p ■ n) 2 /n 2 , 
implying that the differential operator d/dp + can be replaced by d/dn a using a chain rule, 

d n d 

P + TT ( f )= Va T~ 0> ( 7 ) 

dp + v ■ n dn a 

where v = (1,0,0) is a dimensionless vector along p. The operator d/dn a applies to , 
leading to 



t) 2 A 

v a —N^ = v a {N^ a l v + N QU F) , (8) 

v ■ n dn a 

with the special vertex 

V a = 7 • 9 

v ■ nn ■ I 

The momentum Z M (l u ) is contracted with a vertex the differentiated gluon attaches, 
which is then replaced by a special vertex. For each type of vertices, there exists a Ward 
identity, which relates a diagram with the contraction of Z M (l u ) to the difference of two 



diagrams [HJ. A pair cancellation then occurs between the contractions with two adjacent 
vertices. Summing diagrams with different differentiated gluons, the special vertex moves 
to the outer end of a parton line. We arrive at the derivative, 

described by Fig. 1(a), where the square in the new function represents the special vertex 
v a . The coefficient 2 comes from the equality of the two new functions with the special 
vertex on either side of the final-state cut. 

To obtain a differential equation of <ft, we need to factorize subdiagrams containing the 
special vertex out of 0. The factorization holds in the leading regions of the loop momentum 
I that flows through the special vertex. The collinear region of / is not leading because of 
the factor l/(n ■ I) in v a with nonvanishing n 2 . Therefore, the leading regions of I are soft 
and hard, in which the subdiagrams are factorized from into a soft function K and a hard 
function G, respectively. The remaining part is the original distribution function <fi. That 
is, is expressed as the convolution of the functions K and G with 0. 

The lowest-order contribution to K from Fig. 1(b) is written as 

s (x, k T ,p + ) = 4>sv(x, k T ,p + ) + 4> sr (x, k T ,p + ) , (11) 



with 



r d 4 ~ e l v@v v 1 



(x,k T ,p + ), (12) 



(27T) 

'9 2 C F p e I -^iV^(/)^27rz5(/ 2 )0(x, |k T + l T |,p+) , (13) 
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corresponding to virtual and real gluon emissions, respectively. CV = 4/3 is a color factor, 
and 5K an additive counterterm. The ultraviolet pole in Eq. (|T2| ) is isolated using the 
dimensional regularization. To work out the loop integral in Eq. (pT3|) explicitly, we employ 
the Fourier transformation from kx space to b space. The convolution of the subdiagram 
with in the loop momentum Ij- is then simplified into a product. This is the reason the 
kj< resummation should be performed in the space conjugate to k F - The combination of 
Eqs. (|12|) and (|i~3"D then gives 



s (x,6,p+) =K(l/(b f i),a s ( f i))4 > (x,b,p + ) , (14) 

with 

d 4 ~ £ l „ ,„vV 



K = zg 2 C F ^J *—N vfl {iy 



i + 2m5(l 2 )e llT - h 
r 



5K , (15) 



where the factor exp(«Lr • b) is introduced by the Fourier transformation. 

Note that we have applied the soft approximation to the real gluon emission: 

<f>{x + l + /p + , |k T + l T \,p + ) « <f>(x, |k T + l T |,p + ) , (16) 

when writting Eq. (^3|) . This approximation implies that 1/6 serves as an infrared cutoff for 
the evaluation of the subdiagrams shown in Eq. dl5|) . Combined with the hard function G, 



which is characterized by the scale p + } the logarithms \n(p + b) are formed. Hence, Eq. (16) is 
associated with the k^ resummation. We shall show in the next section that the alternative 
soft approximation 

<f>(x + l + /p + , \k T + l T \,P + ) ~ <t>{% + l + /p + , k T ,P + ) , (17) 

implies the infrared cutoff p + /N in moment space. Combined with G, the logarithms ln(l/iV) 
are produced, so that Eq. QT7| ) is associated with the threshold resummation. 
The lowest-order contribution to G from Fig. 1(c) is given by 

(fi h (x, b, p + ) = G(xp + / fi, a s (/i))0(x, b, p + ) , (18) 

in b space, with 



G = -ig 2 C F tf 



r7 + 



(19) 



(xp — I) 2 v ■ I 

where 5G is an additive counterterm. The second term in the above expression, whose sign 
is opposite to that of (f> sv , is a soft subtraction. This term avoids double counting, and 
ensures a hard momentum flow in G. At intermediate x, G is characterized by a large scale 
p + as stated above. 

Using the definition = (f> s + (f> h , Eq. fllCf ) becomes 

d 



K(l/(bfi),a a (ji)) + G{xp + / n,a s (n)) 4>(x,b,p 



(20) 



where the sum K + G is the coefficient function C mentioned before. A straightforward 
calculation leads Eqs. ([15]) and (|19| ) to 



K(l/(6 /U ),« s ( /U )) = ^C F ln-i, 
G(xp + 1 /x, a s (/j)) = ^^CVln 2 ^ 



7T 



(21) 
(22) 



with z/ 2 = (f • n) 2 /|n 2 | being a gauge factor, where we have assumed n 2 < f2[. Constants 
of order unity, which are irrelevant to our discussion, have been neglected. Equation ( |22|) 
confirms our argument that <fi depends on p + via the ratio (p ■ n) 2 /n 2 . For a detailed 
derivation of Eqs. ( pl|) and (|22]) , refer to 0. 

The functions and G possess ultraviolet divergences individually as indicated by their 
counterterms. These divergences, both from the virtual gluon contribution (f) sv , cancel each 
other, such that the sum K + G is RG invariant. The single logarithms ln(6/i) and ln(p + //i), 
contained in K and G, respectively, are organized by the RG equations 



d k d r 

an djj, 



(23) 



The anomalous dimension of K, X K = fidSK/d/i, is given, up to two loops, by [TT 



7* = (-) C F 

7T \ 7T 



C A f5I - t\ - A n/ 
V 36 12 ) 18 ; 



(24) 



with rif being the number of quark flavors, and Ca = 3 a color factor. 



Solving Eq. (p3|), we have 

K(l/ (b/i), a s (fi)) + G(xp + /fi, a s (/i)) 

= K(l, a s (l/b)) + G(l, a s {xp + )) - f" + ^ 7 *(a.00) , 

Jl/b /i 

= - r + -lK(a s ^)) , (25) 

where the initial condition K(l,a s (l/b)) vanishes, and Inu in G(l, a s (xp + )) has been 
dropped, since the gauge factor will be cancelled as computing physical quantities. Inserting 
the above expression into Eq. (|20|) , we obtain the solution 

(j)(x,b,p + ) = A k (b,xp + )(f>(x) (26) 

with the Sudakov exponential from the kx resummation, 

Afc(6, xp + ) = exp 

We have set the upper bound of the variable p to xp + , which corresponds to the final 
condition (f>(x,b,p + ), and the lower bound to 1/6, such that the initial condition <p(x) does 
not contain the logarithms ln(p + b). This statement will become essential, when Eq. (|27|) is 
compared with the exponential from the threshold resummation for 0. 



x p + dp fP d/j, 2 . . . . 

1/6 p Jl/b 2 /I 2 



(27) 



III. THRESHOLD RESUMMATION 

In this section we derive, using the CS formalsim, the threshold resummation for the 
quark distribution funciton, 

<K*,P + ) = I ^-e-^ + y-(p\q(y-)U + q(0)\p) , (28) 

J ATT I 

which is obtained by integrating Eq. (|5|) over kj>. Comparing each step of the derivation, the 
relation of the threshold resummation to the fcr resummation will be clear. The argument 
p + represents the logarithms ln(l — x)p + . According to the same reasoning, depends on 
p + through the ratio (p ■ n) 2 /n 2 , and the chain rule in Eq. (|7|) holds. Following Eqs. @ 
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and (|T0|) , we have the factorization of into the convolution of subdiagrams containing the 
special vertex with 0. Similarly, subdiagrams with soft and hard loop momentum flows are 
absorbed into the functions K and G, respectively. 

The lowest-order contribution to K from Fig. 1(b) is written as 



(p s (x,p + ) = (f) sv (x,p + ) + (j) S r{x,p J ' 



(29) 



with 



ig 2 C F n 
ig 2 C F jjf 



d A ~ e l 
(27r) 4 - £ 

d 4 ~n 



DV 1 



5K 



{x,p + ) 



jc^N u p{l)^2m5(l 2 )<P{x + l + /p + ,p + ) 



(30) 
(31) 



corresponding to virtual and real gluon emissions, respectively. The virtual gluon contribu- 
tion is the same as that in Eq. (0), and Eq. (0) is the consequence of the soft approximation 
in Eq. ([l?]). If retaining transverse degrees of freedom of a parton at the beginning, Eq. (|17D 
allows us to integrate out the k T dependences of 4> sr and of </>, and we still arrive at Eq. ( j3T|) . 
Inserting the identities 

f 1 d£5(£ -x) = l, t di8{i -x- l + /p + ) = 1 , 

J X J X 

into Eqs. (|30|) and ([H]), respectively, Eq. (|2~9"D is reexpressed as 



(32) 



(33) 



with 



d 4 ~ £ l 



(2n 



K = ig 2 C F ^ / 1 ^-N vp {l) 
+2m5(l 2 )5 ( 1 - ~ - 



6(1 -x/t) 



P 



x 



5K5 1 - - 



To obtain the above expression, we have adopted the approximation 



(34) 



5 1 



x 



which is appropriate in the considered region with x 



(35) 
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To work out the £ integration explicitly, we employ a Mellin transformation from mo- 
mentum fraction (x) space to moment (N) space, 

s (iV,p+) EE ^ dxX N - 1 s (x,p + ) , 

Jo 

= K(p + /(N f i),a s ^)) ( p(N,p + ) . (36) 

with 

K(p + /(Nfi),a s (fi)) = f dzz N - 1 K{{\ - z)p + /fi,a s (fi)) . (37) 

J 

The convolution between K and <ft is then simplified into a product. This is the reason the 
threshold resummation should be performed in moment space. The Mellin transformation 
for the threshold resummation and the Fourier transformation for the k F resummation then 
play the same role. As shown in the Appendix, K is given by 

K(p + /(Nf,), a s (»)) = ^C F ( f 1 dz Z ^^ + ln^l , 

7T \J0 1 — Z /i J 

= ^tl CFi np, (38) 

with the same counterterm 5K. To derive the second expression, we have identified the 
integral over z as ln(l/iV), which is valid up to corrections suppressed by 1/N. Hence, soft 
real gluon emissions produce the logarithms ln(p + /iV) as stated before. 

The lowest-order contribution to the hard function G from Fig. 1(c) is written as 

MN,P + ) = G(p + / M ,« s (/x))0(iV,p + ) , (39) 



in TV space, where the expressions of G have been given in Eqs. ( |19D and (|22|) with x = 1. In 
conclusion, the functional forms of K and G in the threshold resummation are the same as 
Eqs. (|2ll) and (|22|) , but with the scales 1/6 replaced by p + u/N and xp + by p + , respectively. 
Using = (p s + (j>h, we have the differential equation 

P + ^(N,P + ) = 2 [K(p + /(N^,a s (fx)) + G(p + /fi,a s ^))} 0(iV,p + ) . (40) 

Again, K and G contain ultraviolet divergences individually, but their sum K + G is RG 
invariant. The RG solution of K + G is given by 
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pp + du 

Jp+/N [I 

To sum ln(l/iV) by means of Eq. (0), we make the replacement 

+ d(f) p + d(f) 
P dp+ = ~Nd(p+/N) ' 

That is, the characteristic scale p + of G is frozen, when solving Eq. 



(41) 



(44) 



(42) 

We then obtain 

<P(N,p + ) = A t (N,p + )<p^, (43) 

with the exponential from the threshold resummation, 

A t (N,p + ) = exp - / -f JL. 7K ( a .Qi)) 
jp+ p jp^ /j, 

which exhibits enhancement. 

We have set the upper bound of the variable p to p + /N, which corresponds to the 

final condition <ft(N,p + ) with p + /p = N, and the lower bound to p + , such that the initial 

condition cj)^ does not contain ln(l/iV) because of p + /p = 1. In this way the N dependence 

of 4> is grouped into the exponential A t . Contrary to Eq. (|27D, the lower bound of p is larger 

than the upper bound in At. Hence, the kx resummation gives the evolution of a parton 

distribution function from the small scale 1/6 to the large scale p + , while the threshold 

resummation gives the evolution from the large p + to the small p + /N. Since the directions 

of the double-logarithm evolutions are opposite, their resummation effects are also opposite. 

Employing the variable change p — (1 — z)p + and pi = \^Xp + , Eq. d44]) becomes 

•i-i/tf dz r 1 dX 



A t (N,p + ) =exp 



-j K (a s (\/Xp + )) 



1 — Z J(l-z) 2 X 

It can be easily justified that the above expression is equivalent to 



(45) 



A t (N,p+) =exp 



r N-l 



dz- 



l r 1 



dX 



/() 1 — 2 .7(1-2)2 A 

up to 0(1/N) corrections, which has appeared in ||. 



1k{u s {^Xv + )) 



(46) 



IV. THRESHOLD RESUMMATION FOR A JET 
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As stated in the Introduction, the scattered quark in DIS carries the momentum p' = 
C,p + q = ((£ — x)p + ,q~,0), where ^p is the initial quark momentum and q the photon 
momentum. In the threshold region with ar — »• ^ — > 1, p 7 possesses a large minus component 
p'~ = q~ but a small invariant p' 2 = 2(£ — x)p •# = (£ — This scattered quark 

produces a jet of particles in this region, to which involved radiative corrections contain 
both collinear divergences from I parallel to p' and soft divergences from I — > 0. Because of 
the cancellation of soft divergences, these corrections are mainly collinear, and absorbed into 
a jet function J. This jet function J depends on the ratio p' 2 / s = £ — x, and more precisely, 
on the momentum fraction w = 1 — x/£ in a factorization formula. For convenience, we 
parametrize the scattered quark momentum as p' = (wp'~ ,p'~ , 0), which can be achieved by 
choosing a frame with p + = q~ = p'~. 

We perform the threshold resummation for J(w,p'~) in the CS framework. It will be 
confirmed below that the large scale p'~ appears in the ratio 

(n ■ v') 2 



n 2 p' 2 



(47) 



in the case with nonvanishing jet invariant mass, where v' = (w, 1, 0) is a dimensionless vector 
along p' . The above ratio differs from (p-n) 2 /n 2 in the massless case for the quark distribution 
function, which has been discussed in the previous section. We shall demonstrate that this 
difference turns the threshold resummation from enhancement for a parton distribution 
function to suppression for a jet function. Note that Eq. (f47|) does not approach (p' ■ n) 2 /n 2 
as w — > 0. This explains why the threshold resummations for a jet function and for a parton 
distribution function do not coincide with each other in the large N (small w and large x) 
limit. 

The chain rule for J, corresponding to Eq. (|7|) for <f), becomes 

- p >-^—j = — v 'J—j. (48) 

dp'~ v ' ■ n a dn a 

Obviously, the extra minus sign on the left-hand side of the above equation will lead to 
suppression. Following the same procedures as in Sec. Ill, we obtain the derivative 
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-p'--^-J(w : p'-) = 2J(w,p'-) , (49) 



where the new function J contains the special vertex 



2 / 

n v ' 



V'a = ~ —i ■ (50) 

v' ■ nn ■ I 

The coefficient 2, again, comes from the equality of the two new functions with the special 
vertex on either side of the final-state cut. 

Similarly, the leading regions of the loop momentum I flowing through the special vertex 
are soft and hard, in which subdiagrams containing the special vertex are factorized out of 
J into a soft function K and a hard function G, respectively. The lowest-order subdiagrams 
for K and for G are basically the same as those in Figs. 1(b) and 1(c). The contribution 
from Fig. 1(b) is written as 

J s {w,p'~) = J sv (w,p'~) + J sr (w,p'~) , (51) 



with 



r d 4 ~ e l v lf3 v' u 1 



J(w,p'-) , (52) 



J sr = ig 2 C F J ^^N uP {l) V -^2^5{l 2 )J{{p' - l) 2 /s,p>-) , (53) 

corresponding to virtual and real gluon emissions, respectively. The virtual gluon contribu- 
tion is the same as in Eq. ([!(]) except for the substitution of the vector v' for v, and thus the 
counterterm is also SK. The first argument of J associated with the real gluon emission, 

(p> -If /+ wl- 

~ = w ~ — ~ — , 54 

s p p 

differs from the corresponding one in Eq. (fJTp. 
Inserting the identities 

r w r w ( l + wl~\ 

J d V^ w -y) = 1 ' J dy8 \w - y - — - —j = 1 , (55) 

into Eqs. (|52|) and (^), respectively, Eq . (|5lD is reexpressed as 

14 



J s (w,p'-) = H-^-K (^-JL^ aa ^)) J(y, P >-) , (56) 
Jo l-y \l-y fj, J 



with 

d 4 ~ € l ,„v^v v 



+2m5{l 2 )5 



1 (w-y 



2 . r , w — y l + wl 



1 — y p'~ p'~ 



To obtain the above expression, we have adopted the approximation 

(w — y l + wl~ \^ (w — y l + wl~\ . . 

which is appropriate in the considered region with w — > 0. 

To work out the y integration explicitly, we apply a Mellin transformation from momen- 
tum fraction (w) space to moment (N) space, 

J s {N,p'-)= f 1 dw{l-w) N - 1 J s {w,jf-), 
Jo 

= f-^- f d W (l- W ) N -iK(^ P —,a s (»)) J(y,p'-), 
Jo l-y J y \ 1 - y /i J 

(59) 

Using the variable change w = z(l — y) + y, the above expression reduces to 

J s (N,p'-) = K{j/-/(N^,a a ^))J(N,p'-) , (60) 

with 

K(j/-/(N»), = f 1 dz(l - z) N - l K(zp'~/^ a s (/i)) . (61) 

Jo 

Replacing the variable z by 1 — z, and following the procedures in the Appendix, the result 
of K is 

K(p'-/(Nt,), = ^±C F ( f 1 dz zN ~ l ~ l + In ^) , 

7T \Jo 1 — Z V fj, J 

= ^C Fi nfr, (62) 

71 Nl/'fJ, 
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with the gauge factor v' 2 = (n ■ v') 2 /\n 2 \, which confirms that J depends on p'~~ via the ratio 
(n ■ v') 2 /(n 2 p'~ 2 ) in the massive case. 

An alternative Mellin transformation adopted in || is 



oo 

J s (N,p'~) = I dwe- Nw J s (w }P '-) , 
Jo 

poo roo 

= dy dwe- Nw K((w-y)p'-/^a s (fi))J(y,p'-). (63) 

JO Jv 



Employing the variable change w = z + y, the above expression reduces to Eq. (|60|), but 
with 



oo 

K(p'-/(N tJl ),a s (n))= I dze- Nz K(zp'-/^a s (fi)) . (64) 



o 



It is easy to justify the equivalence between the transformation in Eq. (|63"D and that in 
Eq. ©. 

The lowest-order contribution to the hard function G from Fig. 1(c) is given by 

J h (w,p'~) = G(y/wp'~ /p,,a s (p,))J(w,p'~) , (65) 
with the hard function 



r r/ 4_e / v'P 



y-/ v 2 P " 

r7 - 



5G . (66) 



(p'-lf p v-2p'-l_ 

Note the modified eikonal approximation for the quark propagator (p' — I) 2 w p' 2 — 2p' ■ I, 
since we are considering the massive case with p' 2 ^ 0. As shown in the Appendix, the 
expression for G is 

G = -^ f ln^. (67) 

As applying the Mellin transformation in Eq. ( |63|) to both sides of Eq. (p5"l), we simply 
replace the momentum fraction w in G by 1/N because of 

roc _ 1 roo _ 

/ dwe- Nw lnwJ(w,p'-) =ln— dwe~ Nw J(w,p'-) , (68) 
Jo N Jo 

which is valid up to 0(1/ N) corrections. Equation (|B"5|) then leads to 
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J h {N,p'-) =G(p'-/(VNfi), a a (jj,)) J(N,p'-) . 



(69) 



Using J = J s + Jh, Eq. (f49|) becomes 

-p'-^A^p'-) = 2[K(p'-/(N^,a s (ix)) + G(p'-/(VN^,a s (fi))]J(N,p'-) . (70) 
The RG analysis gives 

, 1 rp'- 2 /N An 2 

K(p>-/(Nri,a s (ri)+G(p'-/(VN»),a s (ri) = -- a/ a -^J K (a s ^)) . (71) 

Solving the differential equation ( [70|) under a modification of the derivative similar to 
Eq. ([42]), we derive 



(72) 



7K-(a a (/i)) 



(73) 



J(iV,p'-) = A;(iV ! p'-)j( ), 

with the exponential from the threshold resummation, 

A t (N,p ) = exp / / 

We have set the upper bound of the variable p to p'~ /N and the lower bound to p'~, which 
correspond to the final condition J(N,p'~) and to the initial condition respectively. 

At last, adopting the variable changes p = (1 — z)p'~ and fi = \/Xp'~ , Eq. (|73f ) is 
reexpressed as 

v (/: r 1 -* r/A 

l(l-z) 2 



A' t (N,p'-)=exp 
which is equivalent to || 



/ ; / -T-lK\pt s \y\p )) 

J0 I - Z J(l-z) 2 X 



(74) 



A[(N,p')=e W 



'(l-z) 2 

up to corrections suppressed by 1/JV. 



Jo 1 — z Jn-z) 2 X 



(75) 



V. UNIFIED RESUMMATION 



In the previous sections we have shown that both the kx and threshold resummations can 
be reproduced in the CS framework by introducing the transverse and longitudinal cutoffs 
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for soft real gluon emissions, respectively. It is then possible to derive a unification of the 
two resummations for a parton distribution function in the same formalism by keeping both 
cutoffs. Consider the quark distribution function defined by Eq. (^). Adopting the similar 
reasoning, we have Eq. (|T0D, and the factorization of <fi into the convolution of subdiagrams 
containing the special vertex with 0: the subdiagrams with soft and hard loop momentum 
flows are absorbed into the functions K and G, respectively. 

The lowest-order contribution to <p s from Fig. 1(b) is written as Eq. (|ll|), with 
4>sv(x, kT,p + ) the same as Eq. (|T2|) and 

sr = ig 2 C F ^ ( 7 ^LN uP (l)^27ri5(l 2 )<j ) (x + l + /p + , |k T + l T \,p + ) . (76) 

Note that both the l + and It dependences of <fi in the integrand of <p sr have been retained. 
To work out the It integration explicitly, we employ the Fourier transformation from fcy 
space to 6 space as shown in Eq. (|T5|). Following the derivation of Eq. (|34]), <p s becomes 

<M*,6,p+) = f^K ((l- |) «.(/*)) 0M,P + ) , (77) 



with 

~8(l-x/t) 



+M5{P)8 (l - | - l l\ e^ b 



I 2 



SKsll--) . (78) 



Further applying the Mellin transformation to Eq. (77), we obtain 

UNAP + ) = K(p + /(N t i),l/(b»),a s (ti))<f>(N,b,p+) , (79) 

with 

K(p+/{Nri, 1/(6//), a,(fi)) = C dzz^Kdl - z)p + /», 1/(6//), a a (ji)). (80) 

Jo 

The convolutions between K and in l + and It are then completely simplified into a 
multiplication under the Mellin and Fourier transformations, respectively. 
It is easy to work out the integration in Eq. (j30|), and the result is 



K(p + /(Nn),l/(bri,a.Qi)) = 



71 



6// °\ N / 



(81) 



Kq being the modified Bessel function. We examine the large p + b and N limits of the above 
expression. For p + b ^> N, we have K$ — > and 

1 



(82) 



which is the soft function for the resummation in Eq. fl2T|) . For N p + b, we have 
K « — \n(up + b/N) and 



^ -» — C F ln-— 

7T iV// 



(83) 



which is the soft function for the threshold resummation in Eq. (|3^). Hence, Eq. (]81f) is 
indeed appropriate for the unification of the k F and threshold resummations. 

The lowest-order contribution to G from Fig. 1(c) is the same as in Eq. fl3~9"|). Using 
4> = 4> s + 0ft, Eq. ( |I0|) becomes 

d 



v 



dp^ 



4>(N,b,p+) = 2 K(p+/(Nfi), 1/(6//), a s (//)) + G(p+//W/*)) <!>(N,b,P + ) - (84) 



As solving the RG equations (§3), we allow the variable // evolves from the characteristic 
scale of K to the scale of G. Equation (^) implies the characteristic scale of order 



iexp 



' p + b^ 



(85) 



for the unified resummation. We discuss the cases for p + b ^> iV and for N ^> p + b first, 
which will help the derivation of the unified resummation. The solution of K + G is written 

as 



K(p + /{Nfi), 1/(6//), a s {fi)) + G(p+/fji, a a (p)) 
- [ P — IkM/j)) , for p + b^>N , 

71/6 // 

-^7A-(a s (//)) , for iV > p+6 , 

P+/AT // 



(86) 
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indicating that the distribution function <p(N,b,p + ) involves ln(p + b) and ln(l/iV) in the 
p + b — > oo and iV — > oo limits, respectively. This can be easily understood by ignoring the 
variation of 7#, and performing the /i integration directly. 
Inserting Eq. (^) into fl34|), we obtain the solutions 

dp fP d\x 



cj)(N, b,p^ 
<f>(N,b,p^ 



exp 
exp 



2 / — / — IkMii)) 

Jexp[-K (p+b/N)]/b p Jl/b fi 
fP + dp fP dfX ' 

-2 / — / — 7M<w)) 

Jexp[-K (p+b/N)]/b p Jp+ fl 







(87) 
(88) 



for p + b ^> N and 3> p + b, respectively. We have chosen the characteristic scale in Eq. fl85|) 
as the lower bound of the variable p. To unify the above expressions, we replace the lower 
bounds of \x by 



^ex V [-K (p + b)] 



(89) 



that is motivated by Eq. (|85|) . At last, the unified resummation is given by 

0(iV, b,p + )=A u (N,b,p + )^°\ 



(90) 



with the exponential 
A u (N,b,p + ) = 



exp 



-2 / — / — 7M<W)) 

Jexp[-K (p+b/N)]/b p Jexp[-K (p+b)]/b jJ, 



(91) 



which is appropriate for arbitrary p + b and N. It is easy to justify that Eq. ( |9~T{ ) approaches 
the fcy resummation in Eq. (E7p as 6 — > oo, and approaches the threshold resummation in 
Eq. (H) as b -> 0. 



VI. CONCLUSION 

In this paper we have explored the relation between the and threshold resummations, 
and discussed the difference between their applications to a parton distribution function 
and to a jet function in the CS framework. It has been understood that the summations 
of the logarithms \n 2 (p + b) and ln 2 (l/7V) are determined by the soft approximations for real 

20 



gluon emission in Eqs. ( |16D and (|17D, respectively. Suppression from the kx resummation, 
as indicated by A& in Eq. (|27j), and enhancement from the threshold resummation, as indi- 
cated by At in Eq. (pjQ, for a parton distribution function is the consequence of the opposite 
directions of their double-logarithm evolutions. While the switch of the threshold resumma- 
tion from enhancement for a parton distribution function to suppression for a jet function 
is attributed to the change from a massless case to a massive case. The kr resummation 
for a jet function remains suppressive as shown in Eq. @), since the jet momentum and the 
parton momentum are in the same form. Hence, we do not repeat the derivation here. 

We have also unified the kx and threshold resummations by keeping simultaneously 
the longitudianl and transverse loop momenta in the parton distribution function 0(x + 
l + /p + , |k T + lx\,p + ) as shown in Eq. ([761). Comparing the result in Eq. (|91|) with the kx 
resummation in Eq. (f2~7|), which gives suppression, and with the threshold resummation in 
Eq. (|44|), which gives enhancement, the unified resummation exhibits both behaviors: it is 
suppression in the large b region (p + b ^> N), and turns into enhancement in the small b 
region (N ^> p + b). That is, Eq. (^TJ) displays the opposite effects of the kx and threshold 
resummations at different b. The behavior of the unified resummation can be understood 
as follows. For an intermediate x, virtual and real soft gluon corrections cancel exactly in 
the small b region, since they have almost equal phase space. Hence, there are only single 
collinear logarithms, namely, no double logarithms. In this case the Sudakov exponential 



approaches unity as b < l/p + fL2l , indicating the soft cancellation stated above. However, 
at threshold (x — > 1), real gluon emissions still do not have sufficient phase space even 
as b — > 0, and soft virtual corrections are not cancelled exactly. In this case the double 
logarithms ln 2 (l/iV) persist and become dominant. Sudakov suppression then transits into 
enhancement, instead of unity, as b decreases. 

The unified resummation obtained in this work has important applications to studies of 
some QCD processes, such as dijet production In this experiment one jet (the trigger 
jet) is required to be in the central rapidity region, while the other jet (the probe jet) has 
any rapidity up to 3.0. With the large rapidity, dynamics of a hadron at higher x values is 
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probed, so that the threshold resummation is necessary. On the other hand, the differential 
dijet cross section versus the transverse energy of the trigger jet is measured, which further 
demands the inclusion of the resummation. This subject will be discussed elsewhere. 

This work was supported by the National Science Council of Republic of China under 
Grant No. NSC-88-2112-M-006-013. 



APPENDIX 

In this Appendix we present the details of the loop calculations involved in this paper. 
(1) Equation (|38|) 

The virtual gluon contribution to the soft function K associated with the threshold 
resummation for the quark distribution function is given by 

K. = M J - SK , (92) 

where the infrared regulator a is introduced for convenience, and will approach zero at last . 
The first term g v p in N v p gives a vanishing contribution because of v 2 = 0. The contribution 
from the second term —n v lp/n ■ I cancels that from the fourth term n 2 l v lp/(n ■ I) 2 . Hence, 
we concentrate only on the third term —npl v /n ■ I, which leads to 

The above integral has been evaluated in ||, and thus we quote the result directly: 

K, = ?Mc F \n a -, (94) 

7T \i 

with the counterterm 5K = —asCp/ijie). 
The soft real gluon emission contributes 

= 3 * cf i iSf^ 2is{ ' 2 - a2)s{i - : - ' +/p+) ' (95) 

22 



with z = x/£, where e has been set to zero, since the integral is ultraviolet finite. Performing 
the integrations over l~ and l F , Eq. (p5|) reduces to 



7r (1 — — a 2 /(2p + is) 2 
with z/ 2 = (f • n) 2 /\n 2 \. Combined with Eq. (@), the function K in moment space is, 
according to Eq. fl37D, written as 



K( P y(Nri,p+) = ^C F 

7T 



i , z N - 1 (l-z) a 
dz h In — 

(1 — z) 2 — a 2 /{p + v) 2 \i 



(97) 



A simple identity z 1 = (z N ~ l — 1) + 1 isolates the logarithmic term In a in the z integral, 
and the above expression becomes Eq. fl38|) . As expected, the infrared regulator a has 
cancelled between the virtual and real gluon contributions. 
(2) Equation (||) 

We need to evaluate only the real gluon part, 

K r = ig 2 C F J ^N uP {l) V -^2Tti5{l 2 - a 2 )5{z - 1+ /p'~ - wl~ /p'~) , (98) 

with z = (w — y)/(l — y). The first term g v p in N v p gives a negligible contribution because 
of v' 2 = 2w — > 0. The contribution from the second term —n v lpjn ■ I cancels that from the 
fourth term n 2 l v lj3/{n ■ I) 2 . Hence, we consider the third term —n^l^/n ■ I, which leads to 

K r = g 2 C F J j0y 3 ^-jyA 12 ~ ^ - 1+ /P" ~ Wl ~/P'-) ■ (") 

Performing the integration over l F first, we obtain 

r dl+HI' n 2 

K r = g 2 c F / t^-JL-s(z - i + Ip'~ - wi~/ P n , (ioo) 

with the constraint 2l + l~ = l T + a 2 > a 2 . Here we have assumed that the gauge vector n 
possesses only the plus and minus components, i.e., n = (n + ,n~.0) for convenience. Using 
the relation l + = zp'~ — wl~ from the 5-function, the constriant determines the bounds of 
l~ in the integral 

K r = ^c F dr— ^- — - , (ioi) 

2tt Jr. n [n+l~ +n-(zp'- - wl-)} 2 v ; 
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with 

zp'~ — V z 2 p'~ 2 — 2wo? 



mm r» . ' 



2w 



zpf + y/z' 2 p'-' 2 - 2wa 2 

'max" 2w ■ UUZj 



The I integration gives 



K = a s {y) c wn 2 



it rv — wn 

x 



(n+ + wn~)z — {n + — wn z 2 — 2wa 2 /p'~ 
1 



(n + + )z + (n + — wn )\J z 2 — 2wa 2 /p' 2 
After a simple algebra, the above expression is simplified into 



(103) 



7T 2r + (n + — ) z a z /{n z p z ) 

Since the infrared regulator a will approach zero £it lctst, we drop a in the numerator, which 
does not affect the infrared structure of K r . According to the definition 



(n + — wn ) 2 (n • v') 2 /2 (n • t>') 2 

"i 2 nr\ 2 />-) 2 



- i/ 2 « ^ -J- = -z/ 2 (105) 



because of n 2 < 0, Eq. (|104|) becomes 

q B (//) z 

K r = C F — . (106) 

7r z z — v u a A jp' z 

Note that we should neglect the variable w in z/ in order to decouple the integrations over 
y and over w completely in Eq. (|59|). 

Combined with the virtual gluon contribution in Eq. (p4|), we derive in moment space, 



K r (p>-/(Nri,p'-) = ^-C F 

7T 



(1 — z) 2 — v' 2 a 2 jp'~ 2 fx 



(107) 



where the variable change 2; — > 1 — z has been inserted. A similar manipulation isolates ln a 
in the z integral, and the above expression leads to Eq. 
(3) Equation (0) 
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We evaluate first the soft subtraction term in Eq. 



r d A ~ e J v'P 9n' 1 ' 



■2^ e n2 f d*-H 2v'-lp'-n ,, no . 

= —iq CpUL / 7 r- A 7 7Xn , 0/ n - 8G . (108) 

y ^ v' -nJ (2vr) 4 - e (n • Z) 2 / 2 ^' 2 - 2p' ■ I) v ' 

We assume that the coefficient n + of Z~ in the denominator (n ■ I) 2 = (n + l~ + n~l + ) 2 is 
negative. Since the coefficient —2wp'~ of l~ in the denominator p' 2 — 2p' -l = p' 2 — 2wp'~l~ — 
2p'~l + is also negative, a pole from the denominator I 2 = 2l + l~ — l\ exists in the V~ plane 
only for l + > 0. Performing the contour integration over l~ around the pole l~ = l 2 -,/(2l + ), 
Eq. ( |ID^ ) gives 

n n m 2 e f d 2 ~H T r°° Z+(2Z+ 2 + wl 2 T ) , inn . 

G s = 2q z C F n z ii e / - — — — / dl + - — , 7 „ Nn/ , , „ — -r- - 5G . (109) 

y P J (27r) 3 - £ io {2n~l+ 2 + n+l 2 T ) 2 {2l+ 2 -2wp'-l+ + wl 2 T ) V ; 

The integration over l + leads to 

n - in t f d2 ~ H T {n + -wn-) 2 

° s ~ 9 CFf * J (2tt)3- Z 2 (n+ - wn-) 2 + nVf ~ * G ' (U0) 
At last, performing the integration over Z-p, we derive 

G s = -^lB^, (111) 

with the counterterm 5G = o> s Cf / (ire). 

Next we evaluate the first term in Eq. (|66D , 



Gh ~ 19 Cf v~u J kkW 7 ? a ' ( } 

Note that we have set e to zero, because the integral is ultraviolet finite. It is easy to observe 
that poles of Z + exist in the denominator I 2 for l~ < and in (p' — I) 2 for l~ < p'~ . Following 
the similar procedures, we can work out Eq. (|112 ), though the calculation is much more 
complicated. The final expression of Gh is 

G h = ^C F lnw, (113) 



where constants of order unity have been dropped. Combining Eqs. ( Jill ) and ( 113[ ), we 
arrive at Eq. (|67|) . 
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A point that needs to be mentioned is the treatment of the gamma matrices appearing 
in the numerator of Eq. (|112|) . We express 7 + and 7" in terms of fa' and fa: 

, n + fa' — wp'~ fa _ p'~ fa — rT fa' 
p'~(n + — wn~) ' p'~{n + — wn~) 

Since the jet function contains the matrix structure proportional to fa', the products 7+ fa 
and 7~ fa are written as 

, 2n + p' ■ n — wp'~n 2 _ p'~n 2 — IrTp' ■ n 

7 fa = -, ; r , 7 fa = -, ; r- , (H5) 

p'~{n + — wn~) p'~{n + — wn~) 

where we have used the relation /p 1 fan — 2p' ■ n— fa /p 1 ', and neglected the second term, 
which leads to a negligible invariant p' 2 = wp'~ 2 ~ 0. 
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Figure Captions 

Fig. 1. (a) The derivative p + d(f>/dp + in the axial gauge, (b) The 0(a s ) function K. 
The 0(a s ) function G. 
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